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ABSTRACT 



Context. The strength and direction of migration of low mass embedded planets depends on the disk's thermodynamic state. It has 
been shown that in active disks, where the internal dissipation is balanced by radiative transport, migration can be directed outwards, 
a process which extends the lifetime of growing embryos. Very important parameters determining the structure of disks, and hence 
the direction of migration, are the viscosity and the adiabatic index. 

Aims. In this paper we investigate the influence of different viscosity prescriptions (a-type and constant) and adiabatic indices on 
disk structures. We then determine how this affects the migration rate of planets embedded in such disks. 

Methods. We perform three-dimensional numerical simulations of accretion disks with embedded planets. We use the explicit/implicit 
hydrodynamical code NIRVANA that includes full tensor viscosity and radiation transport in the flux-limited diffusion approximation, 
as well as a proper equation of state for molecular hydrogen. The migration of embedded 20 Mem\> planets is studied. 
Results. Low-viscosity disks have cooler temperatures and the migration rates of embedded planets tend toward the isothermal limit. 
Hence, in these disks, planets migrate inwards even in the fully radiative case. The effect of outward migration can only be sustained 
if the viscosity in the disk is large. Overall, the differences between the treatments for the equation of state seem to play a more 
important role in disks with higher viscosity. A change in the adiabatic index and in the viscosity changes the zero-torque radius that 
separates inward from outward migration. 

Conclusions. For larger viscosities, temperatures in the disk become higher and the zero-torque radius moves to larger radii, allowing 
outward migration of a 20-MEarth planet to persist over an extended radial range. In combination with large disk masses, this may allow 
for an extended period of the outward migration of growing protoplanetary cores. 

Key words, accretion disks - planet formation - hydrodynamics - radiative transport - planet disk interactions - adiabatic index 



1. Introduction 

Migration of a low-mass planet embedded in a fully radiative 
gaseous disk can be significantly diff erent from migration in an 
isothe r mal or purely adiabatic di s k (IPaardekooper & Mellema 



2006; 



2008 



Baruteau & Masset 200 



IPaardekooper & Papa oizou 
. .Paardekooper & Mellema! l2008t iKlev & Cridal l2008t 
Kiev et alJ|2009tlAvliffe & Batdl2010l) While all authors agree 
that radiation transport can slow the rate of inward migra- 
tion, there is still a lack of consensus whether the direction 
of migration can be outward. Part of this confusion may be 
due to the sensitivity of the directi on and magnitude of mi- 
gration on global disk parameters (IPaardekooper et al.l l2Q10t 



. _ aldisK p_ _ 

iMasset & Casolil I2010I: IPaardekooper et al.l |2011|). including 



e.g., the radial disk temperature gradient ( Avliffe & Batdl201 ll) . 
Different authors also use different viscosities, usu ally either a 
constant viscosity or a IShakura & Sunvaevl (I1973I) a-viscosity, 
with a typical value of a ranging between 10~ 4 and 10 -2 . 
The simulation of an embedded planet in a global magnetoro- 
tationally unstable accretion disk through direct magnetohy- 
drodynamical calculations including radiative transport is still 
computationally too expensive. However, the disk structure of 
magnetorotationally unstable disks including radiative transport 
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has been calcula ted in the local shearing box approximation 
dFlaig et al.ll2012l) . 

An unperturbed, viscous, fully radiative disk will evolve to- 
wards an equilibrium state, where viscous heating is b alanced by 
radiat ion transport and cooling (as described in, e.g., IKlev et all 
l2009h . This equilibrium state is dependent on the disk mass, 
the viscosity, and the adiabatic index of the gas. Variations 
in viscosity change the radial density and temperature pro- 
files, thus changing the profiles of vortensity and entropy. The 
vortensity and entropy gradients across the horseshoe region, in 
turn, determine the magnit ude and sign of the corotation torque 
dBaruteau & Massed I2008I) . If the torques can remain unsatu- 
rated due to the action of viscosity, the total torque can be posi- 
tive, leading to o utward migration. This effect is also possible in 
isothermal disks (IPaardekooper & Papaloizoul2008l) . 

There is an additional complication to understanding the mi- 
gration of low-mass planets in disks. The rotational states of 
molecular hydrogen are only fully accessible at temperatures 
> 300 K dDecampli et alJll978l : lB"olev et alj|2007l) . As a result, 
the adiabatic index of the gas will transition from y = 7/5 to 
5/3 as the temperature in a disk drops with radial distance. The 
region over which this variation occurs is precisely where we 
expect planetary cores and planets to form in the core accretion 
scenario and begin their initial stages of migration. So far, only 
a constant adiabatic index of y — 1 .4 has been explored in our 
previous simulations with radiative transport. 
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All of these effects on low-mass planet migration are rele- 
vant to understanding whether planet traps can exist, i.e., regions 
in the disk where a protoplanet could experience zero torque. 
Protoplanets migrating from a smaller radius outwards or from 
an outer radius inwards would collect at the zero-torque radius, 
creating areas conducive to planetary mergers, possibly leading 
to the growth of large cores. A planet trap could also be formed 
by surface density changes, which could create a n enhanced 
feeding zone for these cores dMorbidell i et al. 2008); however, it 
is yet unclear how realistic surface density changes are in disks. 
In contrast, radiation transport might allow traps to exist in even 
smooth disk structure. As outward migration is dependent on the 
viscosity, disk mass, and adiabatic index, these parameters will 
influence the radius and breadth of the zero-torque region. In 
this work, we present a study on how these effects modify the 
migration properties of embedded planets. 

The paper is organized as follows: In Section [2] we give an 
overview of our numerical methods. We then describe the in- 
fluence of the adiabatic index, viscosity, and the differences be- 
tween fixed ortho-to-para and equilibrium gas mixtures on the 
disk structure in Section [3] These structural changes influence 
the migration rate of an embedded 20ME ar th planet, which is dis- 
cussed in Section [4] In Section [5] the influence of viscosity and 
of the different gas mixtures on the zero-torque radius are inves- 
tigated. We then summarize and conclude in Section|6] 

2. Numerics and setup 

2.1. Setup 

The protoplanetary disk is modeled as a three-dimensional 
(3D), non-self-gravitating gas. Fluid motion is described by the 
Navier-Stokes equations, where the equations are solved numer- 
ically using a spatially second -order finite volume m ethod that 
is based on the code NIRVANA (IZiegler & Yorkelll997l) . The disk 
is heated solely by internal viscous d issipation, and is allowed 
to coo l by flux-limited diffusion (FLD. lLevermore & Pomraningl 
|1981|) . The FLD approximation allows internally produced en- 
ergy to diffuse radiatively through the optically thick regions of 
the disk and into the optically thin regions, where the energy can 
be radiated away by free-streaming. The flux-limiter interpolates 
between the optically thick and thin regimes and ensures that 
energy loss never exceeds the free-streaming limit. In the code, 
radiative trans port is handled implicitly and it uses t he FA RGO 
(lMassetll2000l) extension as described in iKlev et ail <2009l) . A 
more detailed description of the modeling and the numerical 
methodology is provided in our previous pap ers (IKlev et al.l 
l2009UBitsch & Klevll2010l : lBitsch & Kle\Jl201 lh . 

The three-dimensional (r, 9, <f>) computational domain (with 
266 x 32 x 768 active cells) consists of a complete annulus 
of the protoplanetary disk centered on the star, extending from 
'"min = 0.4 to r max = 2.5 in units of tq = aj up = 5.2AU. For simu- 
lations with planets at larger distances from the central star, the 
computational grid is set to have an inner computational bound- 
ary at OArp and an outer boundary at 2.5rp. In the vertical direc- 
tion the annulus extends 7° above the disk's midplane, meaning 
83° < 9 < 90°. As we assume that the disk is symmetric for the 
upper and lower parts, we apply symmetric boundary conditions 
at 9 = 90°. Here 9 denotes the polar angle of the used spherical 
polar coordinate system measured from the polar axis. The cen- 
tral star has one solar mass M* = M , and the total disk mass 
inside [^^,^3x1 is Mdisk = 0.01M Q . The aspect ratio of the 
disk is calculated self-consistently from the equilibrium struc- 
ture, given by viscous internal heating and radiative diffusion. 



This also determines the surface density gradient in the equilib- 
rium state of the disk. To construct the equilibrium state, we first 
calculate axisymmetric 2D models in the r-9 directions. 

The planet is then embedded at rp = 1.0a j UD . For the 
planet's gravity, we use the cubic potential dKlahr & Klevll200"6T: 
IKlev et alll2009h with r sm = 0.5. The planetary pote ntial and its 
influen ce on migration is discussed in great detail in IKlev et all 
(I2009h . As in our previous simulatio ns we use the opacity law 
provided bv lLin & Papaloizoul (Il985l) . 

2.2. 2D axisymmetric models 

The initialization through an axisymmetric 2D phase (in the 
r-9 plane) reduces the required computational effort substan- 
tially. The evolution from the initial isothermal state towards the 
disk equilibrium between viscous heating and radiative trans- 
port/cooling (henceforth the disk equilibrium state) takes about 
100 orbits, if the disk is started in equilibrium for an isothermal 
gas and constant H/r. The surface density or temperature profiles 
of the initial (isothermal) state are unimportant as the equilib- 
rium state of the disk solely depends on the disk mass, viscosity, 
and adiabatic index of the disk. 

After reaching the disk equilibrium state, we extend this 
model to a full 3D simulation by expanding the grid into the 
^-direction and by embedding the planet. From this starting con- 
figuration it is also possible to investigate the vertical structure 
of an unperturbed disk in the equilibrium state. This corresponds 
to a disk with zero total mass flow because the radial boundaries 
do not allow in- or outflow. 



2.3. Adiabatic index 

The gas in an accretion disk is primarily molecular hydrogen, 
which exists as para and ortho-hydrogen for proton spins that are 
antiparallel and parallel, respectively. If enough protonated ions 
(e.g., H3+) are available to exchange proton spins on timescales 
that are shorter than the dynamical time, the para- and ortho- 
hydrogen should be treated as being in statistical equilibrium 
dBolev et alJ2007h . If the dynamical timescale is shorter than the 
equilibrium timescale, a fixed ratio should be used, where a 3 : 1 
ortho-to-para mix is common for many astrophysical systems. A 
different mixture of gas changes the ratio of specific heats, which 
is the adiabatic index y. The adiab atic index is depen dent of the 
temperature of the underlying gas (Boley et al. 2007). and astro- 
physical disks should transition between y = 5/3 to y = 7/5 
near the frost line. This is represented in Fig.Q] where the adia- 
batic index's dependence on temperature is plotted for the equi- 
librium (blue) and the 3 : 1 ortho-to-para mix (red) gas. The 
adiabatic index can span quite a large region in y for low tem- 
peratures. To make things clearer for the reader, we will refer 
to the equilibrium gas mixture as 0:P(equi) and as 0:P(3:1) for 
the 3:1 ortho-to-para mix. This is also indicated in Fig. Q] In 
our typical disk simulations (with y — 1 .4) the temperature at 
r — 1 .0flj U p is about 65K, which is in a region where the adia- 
batic index changes rapidly, for the 0:P(equi) gas and is close to 
the temperature where the rotational states first activate for the 
0:P(3:1) gas. Such variation could affect an embedded planet's 
migration rate, as the planet modifies the temperature structure 
of its surroundings. For the disks and orbital times explored in 
this study, a fixed ortho-to-para ratio is likely the most physical 
treatment for IT, with the ortho-para equilibration timescale be- 
coming comparable to the orbital timescale outside of 40 AU. 
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Fig. 1. Adiabatic index for the 3:1 ortho to para mix [0:P(3:1)] 
and equilibrium gas fO:P(equi)] configuration in dependence of 
the temperature. 

Nevertheless, we explore both the equilibrium and the 3:1 fixed 
ratio cases to show limiting behaviors. 



2.4. Viscosity 

The viscosity in an accretion disk is typically assumed to be ei - 
ther constant or an a-type viscosity (IShakura & Sunvaevll 19731) . 
for which we use here v = ac 2 s /Q, where c s denotes the local 
adiabatic sound speed and Q the local orbital angular velocity. If 
the accretion disk is in a steady state, the viscosity can be deter- 
mined by the observed mass accretion rate M onto the star: 



(1) 



where M is typically 10~ 9 - lO~ 7 M /yr (lHartigan et alJll995t 
iGullbring et alj 1 1998b and S is the surface density. The vis- 
cosity determines the structure of the disk and thus influences 
the migration of embedded planets. Based on models of accre- 
tion disks, the observed mass accretion rates constrain a to be 
roughly between 10~ 4 and 10~ 2 . As will be discussed later in the 
manuscript, we define the half-disk thickness as 



H = 




(2) 



where c lso is the isothermal sound speed, and P and p denote the 
pressure and density in the disk's midplane. 

2.5. Torque calculations 

In previous work, we have discussed the calculation of the torque 
acting on the planet in great detail. Outward migration seems 
onl y to be possible when the planet is on a nearly circular or- 
bit ( Bitsch & K ev 20101) and close to the midplane of the disk 
dBitsch & Klevl 201 lh . Since we want to investigate the occur- 
rence of outward migration and its implications for planet for- 
mation and evolution, we assume that the planets remain in the 
midplane of the disk and are initialized on circular orbits. Only 
the upper half of the disk (83° < 6 < 90°) is needed for the cal- 
culation, as the lower half is directly symmetric to the upper half 
for these conditions. 

The torques acting on a 20ME a nh planet are calculated to de- 
termine the direction of migration. A positive torque indicates 



outward migration, while a negative torque indicates inward mi- 
gration. As the planet is simulated as a point mass, the plane- 
tary potential needs t o be smoothed. W e use the cubic potential 
dKlahr & Klevll 20061 : 1 Kiev et alJl2009l) for our calculations: 



P ~ 



m,,G 
d 



(£) 4 - 2 (^) 3 + 2 £ 



d 



for d < r sm 
for d > r sm 



(3) 



Here nip is the planetary mass, d = |r - rp| denotes the dis- 
tance of the disk element to the planet, and r sm is the smooth- 
ing length of the potential, usually a fraction of the Hill radius, 
Rn = (l/3m p /M,) ' 3 . The construction of the planetary potential 
is in such a way that, for distances larger than r sm , the potential 
matches the correct l/r potential. For d < r sm , the potential is 
smoothed by a cubic polynomial. The parameter r sm is equal to 
Q.5Ru in all our simulations. 

The gravitational torques acting on the planet are calculated 
by summing over all grid cells, which are treated in this case as 
point masses. We also apply a tapering func tion to exclude th e 
inner parts of the Hill sphere of the planet (ICrida et al.ll2008l) . 
This torque-cutoff is necessary to avoid large, probably noisy 
contributions from the inner parts of the Roche lobe and to disre- 
gard material that is possibly gravitationally bound to the planet 
dCrida et al]|2009l) . Here we assume (as in our previous papers) 
a transition radius of 0.8/^h- 



3. The equilibrium disk structure 

As the adiabatic index is connected to the sound speed, a change 
in the adiabatic index can alter the structure of the disk by modi- 
fying the disk viscosity, affecting the disk equilibrium state. The 
effects of these parameters on the disk structure (without an em- 
bedded planet) are discussed in this section. 

3. 1 . Influence of a constant adiabatic index 

For the following simulations, we limit ourselves to constant val- 
ues of y in the range of 1.05 < y < 1.8, which is larger than 
the range of the adiabatic index in Fig. [TJ We include y values 
smaller than 1 .3 to capture the transition to isothermal disks, and 
higher y values to capture additional possible stiffness resulting 
from magnetic fields. In a first set of simulations we study disks 
with a constant kinematic viscosity of v = 10 15 cm 2 /s. For this 
case, we do not expect a change in the disk structure for different 
values of y because the disk equilibrium state is determined by 
a balance between viscous heating and radiative cooling, both 
of which do not depend of y. Indeed, as can be seen in Fig. [2] 
for a wide range of y, the midplane density and temperature and 
the resulting aspect ratio at r — 1.0aj up for a 2D axisymmet- 
ric disk do not depend on y. Only for small values of y is the 
disk structure different. The reason for this change is the onset 
of convection for y<1.2. As the adiabatic index is lowered, the 
gas becomes more compressible owing to the larger degrees of 
freedom that must be excited before the thermodynamic temper- 
ature can be raised. This leads to a higher midplane density, a 
lower midplane temperature, and a lower aspect ratio (see eq. |2ji 
in the disk at r = 1 .0aj up . 

Radial density and temperature gradients result in vortensity 
and entropy (S oc ■£) gradients. These vortensity and entropy 
gradients across the horseshoe region determine the magnitude 
and sign of the corotation torque, which, if positive, can lead 
to a total positive torque, indicating outward migration in fully 
radiative disks. As the density decreases and the temperature 
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Fig. 2. Densities (top), temperatures (middle), and aspect ratios 
(bottom) from the 2D axisymmetric disks with different (con- 
stant) adiabatic indices, taken after the disks have reached their 
equilibrium state. Each quantity is evaluated using midplane val- 
ues of the disk at a radius r = laj up . The kinematic viscosity has 
been set to v = 10 15 cm 2 /s for all models. 



increases with increasing y (for y < 1.2), the strength of the 
gradient in entropy might change and therefore influence plane- 
tary migration. The influence of this effect might be stronger for 
the 0:P(equi) gas disk, as the adiabatic index can become lower 
compared to the 0:P(3: 1) gas disk. 




Fig. 3. Aspect ratio (H/r) of a fully radiative disk with different 
adiabatic indices without an embedded planet. The aspect ratio 
was computed using midplane values of the constant viscosity 
disks (see eq.|2]i. 



The radial variation of the aspect ratio of a fully radiative 
disk with different adiabatic indices is shown in Fig. [3] 

For y = 1.1, the aspect ratio is smaller at r < 1.5flj up and 
the fluctuations in this region are a sign of convection. However, 
such a small adiabatic index is not predicted for accretion disks, 
as seen by the y(T) profile in Fig. Q] unless efficient heating 
and cooling push the disk toward isothermal behavior or the 
gas is undergoing dissociation or ionization. In contrast, the pro- 
files for a higher adiabatic index show an overall drop in as- 
pect ratio with increasing distance to the central star (after a 
small increase for smaller distances to the central star, r < 0.7). 
Such a feature in the disk's structure is the opposite of what 
one would expect. When taking stellar irradiation into account 
an increasing H/r profile for larger radii may be possible, e.g. 
iDullemond & Dominikl d2004l) . However, the disk's opacity can 
greatly influence the H/R profile. 

With increasing distance from the central star, the temper- 
ature is decreasing. For temperatures less than T ~ 155 K, the 
opacity increases with increasing temperature. At higher temper- 
ature, the opacity decreases. So, with increasing distance to the 
central star, the opacity first increases and, when T ~ 155/iT, the 
opacity decreases. This is exactly the point in the disk (r =s 0.7) 
where the aspect ratio H/r begins to decrease. Please note that 
the height H of the disk is still increasing, but at a lower rate, 
when T < 155 K compared to T > 155 K. Hence, the struc- 
ture of the disk is a result of the opacity law that is used in these 
calculations. 

The H/r profiles for the disks with y > 1.4 are identical, as 
the disk equilibrium state is determined by a balance of viscous 
heating and radiative cooling, neither of which depends on y, at 
least for the constant v assumed for these disks. 



3.2. Influence of viscosity 

As can be seen in Fig. [4] where we display the midplane den- 
sity (top) and temperature (bottom) of fully radiative disks in the 
disk equilibrium state, a change in the viscosity affects the disk's 
equilibrium structure, altering the density and temperature pro- 
files of the disk. In this figure we show the densities and tem- 
peratures for a-viscosities 0.001 < a < 0.008 and for a constant 
viscosity of v = 10 15 cm 2 /s, a value used in our previous sim- 
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Fig. 4. Density (top) and temperature (bottom) profiles of the 
fully radiative disks with different viscosities, all taken in the 
disk equilibrium state. Densities and temperatures are measured 
in the midplanes of the disks. 

ulations dKlev et alJl2009t iBitsch & Kle\Jl201Cl : iBitsch & Kle\i 
I2011I) . 

For increasing ^-viscosity, the density in the midplane of 
the disk becomes smaller, but the differences in density become 
smaller for higher viscosities. 

For the temperature profile, the trend is reversed: a higher 
viscosity results in a higher temperature in the midplane of the 
disk. As the viscosity increases, so does the viscous heating. The 
heating process mainly takes place in the midplane of the disk, as 
the viscosity is highest there. Not only is the temperature higher 
in the midplane, but the temperature profiles also change with 
viscosity. For our constant viscosity of v = 10 15 cm 2 /s, the radial 
slope of temperature was oc r~ L8 at r — laj up , while it now de- 
pends on the chosen a- value. The slope of temperature /3j, where 
T oc r ^ T , is displayed in Table[T]for r = laj up . 

The temperature slope has a d irect influence on the mi- 
gratio n of embedded planets, e.g., iPaardekooper et all (1201 Oi 
1201 lb . as the Lindblad torque and the entropy related corotation 
torque are very sensitive to it. A change in /3t therefore influ- 
ences the torque acti ng on an embedded planet quite severely 
dAvliffe & Batdl201lh . 

The aspect ratio H/r of the disk also increases for increasing 
viscosity (see Table [TJ, indicating another factor that influences 
planetary migration in t hese disks. For isot hermal disks, this has 
been know for a while (jTanaka et al. 20021). 



The density, temperature and H/r profiles of our constant 
viscosity simulation match quite well with the a = 0.005 sim- 
ulations, indicating that one could expect a very similar torque 
acting on a planet embedded in these disks. However, as the a- 
viscosity is not constant in r, one might also expect some dif- 
ferences in the spiral wave densities and temperatures in these 
disks. In addition, the area near the planet might be subject to 
some changes for different viscosities. 

3.3. Influence of a varying adiabatic index 

In Fig. [5] the temperature (bottom) and density (top) profiles in 
the midplane of fully radiative disks are displayed. The disks fea- 
ture different a viscosities and two different ortho-para mixtures 
of the gas [0:P(equi) and 0:P(3:1)]. As expected from simula- 
tions with a constant adiabatic index, an increasing a viscosity 
results in lower densities and higher temperatures in the mid- 
plane of the disk. 

The gradients in density and temperature are comparable to 
the gradients for constant adiabatic indices. It also seems that 
the mixture of the gas does not have a huge influence on the 
temperature and density profiles for large viscosities. However, 
for a = 0.001 there are differences in the density and tempera- 
ture profiles between disks that have different ortho-para ratios. 
For r < 1 .0flj U p, the temperature for the disk equilibrium state is 
lower in the 0:P(equi) than in the 0:P(3: 1). As a result, the den- 
sity is higher in the 0:P(equi) disk. As y is largely reduced for 
temperatures around T ~ 50K for ortho-para ratios 0:P(equi), 
this behavior is consistent with the results of the constant y sim- 
ulations. 

In Fig.|6]the aspect ratios H/r for the same disks are shown. 
The profiles are very similar to those in disks with a constant 
adiabatic index and an a viscosity. In the a = 0.001 case, the 
aspect ratio profile is smaller and approximately constant. For 
the more viscous disks, the aspect ratio increases for small r, 
but then decreases with increasing r continuously. A higher vis- 
cosity leads to a higher aspect ratio in the disk. The 0:P(equi) 
shows a slightly smaller aspect ratio for viscosities a > 0.004 
compared to the 0:P(3: 1). The aspect ratio differs by about 10% 
for the high viscosity cases at r — 1 .0aj up . This difference at this 
location can influence the torque acting on embedded planets, as 
seen in isothermal disks. 

A higher disk vi scosity leads to a larg er temperature in the 
midplane. Recently, IBitsch & Kiev! ( 1201 lh have shown that in 
fully radiative disks with viscous heating, a convective region 
arises near the central star, where its radial extent depends on 
disk mass. A higher disk mass resulted in a larger convective 
zone in the disk. The density and temperature structures of the 



a 




H{r= 1.0aj„p) 


0.001 


1.60 


0.0219 


0.002 


1.70 


0.0277 


0.003 


1.75 


0.0324 


0.004 


1.80 


0.0363 


0.005 


1.80 


0.0393 


0.006 


1.85 


0.0419 


0.007 


1.87 


0.0442 


0.008 


1.90 


0.0468 


v = const. 


1.80 


0.0395 



Table 1. The slope of the temperature profile /3j, where T oc r ^ T , 
and the disk thickness H as a function of disk viscosity. The 
quantities are evaluated at disk radius r — 1 .0aj up . 
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Fig. 5. Midplane density (top) and temperature (bottom) profiles 
for the a viscosity disks with different assumptions for the ortho- 
para ratio of the gas. 




Fig. 6. Aspect ratio H/r of the accretion disks for different a 
viscosities and different assumptions for the ortho-para ratio of 
the gas. 



disk with constant viscosity seem to correspond to a disk with 
a ~ 0.005 (see Fig.©. 

We display in Fig. [7] the velocity in the z-direction for the 
a = 0.001 disks with two different ortho-para gas mixtures, for 
the a = 0.004 and 0:P(equi) disk, and for the fixed y — 1 . 1 disk. 
In black we indicate H(r) and in red the t = 1 .0 line in the disk 
(integrated from the top of the grid toward the midplane). We 
note that when taking both sides of the disk into account, convec- 



tive eddies cross the midplane. However. iBitsch & Klevl (1201 ll) 
have shown that the first signs of convection can also be seen 
in simulations where only one half of the disk is simulated. The 
extent to which the convective region inside the disk depends 
on the adiabatic index should be tested in simulations of disk 
instabilities. As we want to focus here on the influence of the 
adiabatic index on planet migration, we do n ot pur sue this topic 
further. We do note that iRuden & Pollack! (1 1 99 lb have shown 
that a lower adiabatic index can increase the susceptibility of the 
disk to convection. 

For the a = 0.001 disk with either the 0:P(3:1) or the 
0:P(equi) gas, variations in the vertical velocity of the gas in the 
disk's upper layers can be seen in Fig. [7] However, these fluctua- 
tions are limited to the low-density material, and do not reach all 
the way down to the high optical depth region. Convection does 
not appear to be present. Fluctuations such as those in the top two 
panels of Fig.|7]only occur for the thinner disks (small H/r) that 
have a very large vertical range in density, which is difficult to 
resolve numerically (in t his case no-density fl oor is used, possi- 
bly increasing the effect). iNelson et a"fl(l2012l) stated that vertical 
instabilities in isothermal discs can arise, if the viscosity is neg- 
ligible. This vertical shear instability can then drive a non neg- 
ligible transport of angular momentum through Reynolds stress 
with a the corresponding viscous alpha value a = 0.001, which 
is the value of viscosity assumed in the simulations shown in 
the top two plots in Fig. [7] indicating that the shown fluctuations 
might be resulting from the vertical shear instability. However, 
detailed studies of the vertical shear instability for the case of 
fully radiative disks is beyond the scope of the presented work. 

Convection, on the other hand, can only be present in the 
optically thick regions of the disk, where radiative cooling is in- 
efficient. For larger viscosities {a = 0.004), the fluctuations in 
the upper layers mostly vanish, except in the low-density region 
at large r, and are totally diminished for a = 0.008 (not shown). 
Large viscosities suppress the growth of ve rtical oscillations, in 
agreement with the vertical shear instability (Nel son et al.l2012h . 

On the other hand, the fixed y = 1.1 disk shows fluctuations 
reaching all the way down to midplane, where the disk is opti- 
cally thick. Please also note here that the absolute values of the 
velocities are larger for the y = 1 . 1 disk compared with the other 
presented simulations. In this case, we have to check if convec- 
tion is really present in the disk, as we already suspected from 
the aspect ratio profile of the disk (see Fig. [3}. 

In order to test whether convection should be present in the 
disk, we compute the adiabatic gradient as follows: 



d\nT 
dlnp 



(A) 



The adiabatic gradient is evaluated at a two-degree angle from 
the midplane, which roughly corresponds to the half opening an- 
gle of the optical thick parts of the disk. 

For > ySadi, convection can be present in the disk, where 
j8 a di is defined as 



adi 



Tlocal 



1 



Tlocal 



(5) 



Here, yi oca i represents the local adiabatic index, which in our 
case, is the local adiabatic index along the two-degree angle of 
the disk. In Fig. [8] we present the the adiabatic gradients of the 
y = 1.1 disk (bottom in Fig. |7j and the a - 0.004 disk with 
0:P(equi) (3rd from top in Fig. |7}. 

In the case of y — 1.1, /? > ;8 a di, indicating that convection 
should be present in the disk, which corresponds to the vertical 
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Fig. 7. Velocities in z-direction for disks with a = 0.001 and 
0:P(3:1) (top), a = 0.001 and 0:P(equi) (2nd from top), a = 
0.004 and 0:P(equi) (3rd from top), and a fixed y = 1.1 with 
constant viscosity (bottom). The black line indicates H{r) in the 
disk and the red line delineates t — 1 .0 (the optical depth inte- 
grated from the top of the grid toward the midplane of the disk). 
Please note the different color scales and that they only show 
a subset of the range of values to better highlight the internal 
structure. 



fluctuations for this simulation in the optically thick regions of 
the disk. The strong fluctuations in /? indicate that the disk is 
highly unstable. Please also note that in the inner regions of the 
disk, the difference between /? and yS a( ji is about a factor of 5. For 
the a - 0.004 disk, only the inner regions feature (3 > f3 d a. 
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Fig. 8. Adiabatic gradients for the y = 1.1 disk with constant vis- 
cosity (top) and for the a = 0.004 disk with 0:P(equi) (bottom). 
The gradients have been evaluated in the optical thick regions of 
the disk at an angle of 2° from the midplane. 



However, fluctuations in the velocity have not been observed 
(Fig-EJ- In addition, the difference between fi and /J a( ji is at maxi- 
mum a factor of 2, which is much lower than in the y = 1.1 case. 
If convection is present in the a = 0.004 disk, then it must be 
much less vigorous than seen in the y — 1 . 1 disk. 



4. Influence on planetary migration 

Because the above results show that the adiabatic index and the 
viscosity affect disk profiles, we also expect that they will change 
the migration rates of embedded planets. Embedding a small 
mass planet in such a disk will, of course, change the profiles 
of the disk. However, for most viscosities and adiabatic indices, 
the changes that the low-mass planet will cause are limited to 
the vicinity of the planet, so that the overall structure of the disk 
remains intact. To determine the direction of migration of pro- 
toplanets in disks, we embed 20ME ar th planets in our disks and 
measure the torque acting on those planets. A positive torque 
represents outward migration, while a negative torque indicates 
inward migration. The torque acting on planets on circular, non- 
inclined orbits is directly proportional to the migration rate. In 
cases where the planet is embedded in the convective region in 
the disk (for y<l.l), the torques have been averaged over 20 
planetary orbits. 
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Fig. 9. The total torque r tot acting on a 20ME a nh planet embedded 
at r — 1 in v = 10 15 cm 2 /s disks for different (constant) adiabatic 
indices. The torques have been measured in the equilibrium state 
of the disk. 
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Fig. 10. Profiles of the radial torque density T{r) acting on 
20ME a rth planets embedded in disks with different adiabatic in- 
dices (during the disk equilibrium state). In contrast to the total 
torque, the radial torque density has not been averaged in time. 



4.1. Constant adiabatic index 

In Fig. [9] we show the total torque r tot acting on planets that 
are on circular orbits and are embedded in disks that have dif- 
ferent adiabatic indices. The viscosity is held constant at v = 
10 15 cm 2 /s. The torque acting on the planet is positive when the 
adiabatic index is larger than 1.12 , while it is negative, indicat- 
ing inward migration, for y < 1.12. The maximum of the torque 
is at y — 1 .2, and the torque decreases for lower and higher y's. 
The most negative torque is at y = 1.075. 

For y > 1.12, changing the adiabatic index influences the 
overall results of planetary migration only by magnitude, not by 
direction. At high adiabatic indices, the torque acting on the em- 
bedded planet is decreased by a factor of 3 to 4 compared to the 
maximum torque at y — 1 .2. The tendency of outward migration, 
however, remains intact. An increase in y leads to an increase in 
the adiabatic sound speed, which in turn decreases the amplitude 
of the total torque. It also reduces the amplitude of the (negative) 
entropy gradient, which then reduces the amplitude of the (pos- 
itive) corotation torque. This can be seen by the decrease of the 
(positive) spike in the torque density distribution (see Fig. [TOt . 
For lower adiabatic indices, the direction of migration is reversed 
as y approaches the value for isothermal disks. However, the 
torque acting on planet in a radiative disk with y = 1.001 does 
not reflect the torque acting on a planet in an isothermal disk, as 
the disk structur e is different dPaardekooper & Mellemal 120081: 
iKlev et ai]l2012h . On the other hand, the torque has strong fluc- 
tuations in time, indicating that the planet has reached the con- 
vective zone of the disk. 

In Fig.[10l the radial torque density T(r), which is the torque 
exerted by an annulus of width dr and related to the total torque 
by 

r« = J ™ T(r)dr, (6) 

is displayed for 20ME a nh planets embedded in fully radiative 
disks with a constant adiabatic index. All displayed torque den- 
sities F(r) (with y > 1.2) feature a spike just below r ~ 1.0aj up 
on top of the underlying Lindblad torque curve, which is the nor- 
mal pattern for a fully radiative disk with a high adiabatic index 
(IKlev et all2009h . For increasing y, the Lindblad torque and the 
spike in the torque density are reduced, indicating a reduction of 



the positive corotation torque. This reduction in torque density 
is also represented by a decrease in the total torque (see Fig . [9j> ■ 
For y = 1.1, no spike in the torque density is visible, and only 
the Lindblad torque remains. This also is reflected in the total 
torque, which is negative (see Pig . [9j> ■ 

Fig. [TTI shows surface density images of our fully radiative 
disk simulations with an embedded 20ME alt h. The adiabatic in- 
dices for these images are y — l . I (top), y = 1 .4 (middle) and 
y = 1.7 (botto m). We note that th e y = 1.4 case is discussed 
in great detail in IKlev et alJ (l2009h . For the y — 1 . 1 simulation, 
strong variations in the density surrounding the planet are visi- 
ble. The variations are due to convection inside the disk. 

The surface density pattern of the y = 1.7 disk, on the other 
hand, shows no sign of convection. In fact it is very similar to the 
y = 1.4 (middle in Fi g . ITTb . The only difference is that the den- 
sity perturbations in front and behind the planet are not as well 
pronounced. Due to a difference in the adiabatic sound speed, 
which is oc Jy, the opening angles of the spiral waves increase 
with increasing y. The density in the spiral waves is reduced for 
increasing y, resulting in a decrease of the Lindblad torque, as 
can be seen in the torque density distribution (Fig. [Tot , and in 
the total torque (Fig. [9). 

4.2. Influence of viscosity 

In Section [3~!2l we discussed the changes of the structure of plan- 
etary disks with different viscosities. Changes in the temperature 
gra dient lead to changes in th e torque acting on embedded plan- 
ets dPaardekooper et aI1l201 ll) . This sensitivity is highlighted in 
Fig. [12] which shows the total torque acting on an embedded 
20MEanh planet in disks with different viscosities. Planets in 
disks with a < 0.002 (not displayed) feel a negative torque, 
indicating inward migration. A decrease of the torque with de- 
creasing viscosity can also be found in isothermal disks (s ee for 
example our simulations in Fig. 1 in iBitsch & Klevll2O10h . The 
decrease in the total torque is a result of the lower viscosity. If 
the viscosity is very low, the corotation torques saturate and out- 
ward mig ration can not be supported any longer (for isothermal 
disks, see Masset 2001). The planet depletes its coorbital region 
of gas (i.e., begins to open a gap, Fig. [T4b . which reduces the 
torque. 
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Fig. 11. Displayed are the surface density maps for planets on 
fixed circular orbits in fully radiative disks with different adia- 
batic indices (top: y = 1.1, middle: y = 1.4, bottom: y = 1.7) 
and a constant viscosity of v = 10 15 cm 2 /s. The disruption in the 
surface density pattern for the y = 1.1 disk is due to convection 
inside the disk. 



For a viscosities with a > 0.003, the torque acting on 
the planet becomes positive, and for even larger viscosities 
(a > 0.004), the magnitude of the torque begins to converge. 
The torque acting on the planet with constant viscosity of v = 
10 15 cm 2 /s is similar to the one with a ~ 0.006. This is consis- 
tent with the similarities between their temperature and density 
profiles (Fig. HJi. 

In Fig. [13] we display the radial torque density F(r) acting 
on a 20ME art h planet in disks with different viscosity. In the a = 
0.001 simulation, the displayed torque density shows a typical 
Lindblad torque curve, but without the spike at corotation (see 
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Fig. 12. Torque acting on 20ME a nh planets embedded in disks 
with different viscosities (y = 1.4). 



below); a spike in the torque curve at corotation is the typical 
patter n of a fully radiative disk with a high viscosity dKlev et al.l 
l2009h . In fact the torque density is very similar to an isothermal 
disk. 

The a = 0.005 and a = 0.006 torque density patterns 
are consistent with the that of the constant viscosity of v = 
10 15 cm 2 /s. The spike in the distribution near r ~ 1.0a j up is 
nearly identical for the constant viscosity and a = 0.006 simu- 
lations, while the a = 0.005 spike is a little bit higher. However, 
there are some differences. For smaller and larger distances from 
the central star (r < 1 .0aj up and r > 1 .0aj llp ) the torque density 
is larger for the constant viscosity case compared to a — 0.006, 
although all simulations follow the same trend. The underlying 
Lindblad torque distribution is also smaller for the a = 0.005 
disk, due to a difference in the overall disk structure (see Fig. |4). 
The reduction of F(r) for larger a is partly a consequence of the 
increased disk temperature. These small differences are the rea- 
son why the total torque of the constant viscosity simulations is 
somewhat larger compared to the a viscosity disks. 

When a = 0.008, the Lindblad torque is much smaller com- 
pared to the other disks, as the temperature is higher. However, 
the torque density shows the usual spike, indicating a positive 
corotation torque, which can lead to outward migration if it over - 
compensates for the negative Lindblad torque dKlev et al.l2009l) . 
One might suspect from the trend of the simulations that even 
higher viscosities could destroy this effect of outward migration. 
Moreover, modifications to the horseshoe dynamics for large 
viscosities could lim it or cutoff the torque acting on the planet 
(lMassetll200lll2002h . 

In Fig. [14] we display the surface densities of several disks 
with an embedded 20ME a nh planet. Different viscosities are 
shown to gain more insight into the origin of the torque density 
profiles. For very low viscosities (a = 0.001), the planet seems 
to open a very small or partial gap inside the disk. In addition, 
the spiral waves of the planet are more dense compared to the 
constant viscosity simulation (middle panel in Fig. ITTb . which 
may be due to the higher density in the midplane of the disk. A 
low viscosity encourages the formation of gaps in the disk, thus 
explaining the torque density in Fig. [13] 

The a = 0.005 surface density profile seems very similar 
to the constant visc osity profile we discussed in great detail 
in iKlev et al] d2009t) . The density increase ahead of the planet 
(<p > 180° and r < 1.0aj llp ) is clearly visible in both cases, thus 
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Fig. 13. Profiles of the radial torque density T(r) acting on 

20MEanh planets embedded in disks with different viscosities 

(in the disk equilibrium state) and a constant adiabatic index 

(7=1-4). 



creating a nearly identical 'spike' in the torque density distribu- 
tion (Fig. [13). However, the density decrease behind the planet 
(0 < 180° and r > 1.0aj up ) is not so clear in the a = 0.005 
simulation compared to the v = 10 15 cm 2 /s simulation (middle 
panel in Fig. ITTb. This directly reflects on the torque density, as 
it is higher for the constant viscosity simulation compared to the 
a = 0.005 simulation at that distance from the central star, thus 
explaining the higher total torque. 

For a = 0.008, the surface density profile shows the same 
structure as for the constant viscosity simulation, but the sur- 
face density is generally reduced. The spiral waves and the vicin- 
ity near the planet show smaller surface densities. The structure 
ahead and behind the planet is also not as distinctive as in the 
constant viscosity (or a = 0.005) simulation. This all leads to a 
smaller curve in the torque density plot (Fig. [13} and to a smaller 
total torque (Fig. |T2j. 

4.3. Varying adiabatic index 

In addition to different viscosities, we now consider models that 
feature a varying adiabatic index. In Fig.[l5]the total torque act- 
ing on 20MEanh planets in disks is displayed for different viscos- 
ity and gas treatments. For low viscosities (a = 0.001) the torque 
is negative (not displayed), indicating inward migration, and for 
higher viscosities, the torque is positive indicating outward mi- 
gration, as expected from simulations with a constant adiabatic 
index. 

For all shown viscosities, the torque acting on the planet is 
positive, indicating outward migration. In the 0:P(equi) disk the 
torque is higher than the torque for a planet in the 0:P(3: 1) disk 
(with the same viscosity). The largest differences between the 
two ortho -para gas states is observed for the disks with a constant 
viscosity. 

The temperature near the planet is in the region of 807^ < 
Tpianet ^ 130-ST, which is where the separation of the adiabatic in- 
dex between 0:P(equi) and 0:P(3:1) first becomes pronounced 
for decreasing temperatures (see Fig. [T|. This difference in tem- 
perature leads to a change in the total torque acting on the planet 
(see Fig. [9j, which can be up to a factor of two in this temper- 
ature region. Therefore the total torque in the 0:P(equi) disk is 
higher than in the 0:P(3:1) disk. This effect is due to a change 




2.35 & 



Fig. 14. Surface density maps for 20ME art h planets on fixed cir- 
cular orbits in fully radiative disks (with y = 1 .4). The viscosities 
in the disks are (from top to bottom): a = 0.001, a = 0.005, and 
a = 0.008. The snapshots are taken in the equilibrium state. 



in the adiabatic index in the region near the planet. The origin of 
this effect is also described in Section |4~T1 



Fig. Q~6] shows the torque density F(r). As all total torques 
are positive (Fig. [TBI , each simulation shows a spike on top of 
the Lindblad torque in the F(r) distribution just below r ~ lap. 
This spike is largest for 0:P(equi) in the constant viscosity sim- 
ulation, which is also the simulation featuring the highest to- 
tal torque. The spikes in the F(r) distributions are smaller for 
0:P(3:1) disks compared to 0:P(equi) disks, which is also re- 
flected in the total torque (Fig. [TBI . However, the Lindblad 
torque seems strongest for the a = 0.004 disk with the 0:P(equi) 
gas. 
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Fig. 15. Torque acting on 20ME ar th planets embedded in disks 
with different viscosities. The graph features a = 0.004, a = 
0.008, and constant viscosity for both gas ortho-para mixtures. 
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Fig. 16. Torque density of 20ME art h planets embedded in disks 
with different viscosity. The graph features a = 0.004, a = 
0.008, and a constant viscosity for both gas ortho-para mixtures. 



The torque density in the 0:P(equi) disk with a constant vis- 
cosity is slightly smaller around r ~ 1 .0flj U p compared with the 
constant viscosity and y case (see Fig. [13), which leads to a 
smaller total torque. The torque density spike in the disk with 
0:P(3:1) is much smaller than in the constant y case (both for 
constant viscosity again), which results in a smaller total torque. 
The same trend of the behavior of the torque density seems 
to be true for a viscosities with a = 0.008. The spike in the 
torque density is larger for the 0:P(equi) simulations than in the 
0:P(3:1) simulations. This indicates a larger positive corotation 
torque, which results in a larger torque acting on the planet (see 
Fig-OS. This is a consequence of the temperature dependence of 
the adiabatic index, where y has smaller values in the 0:P(equi) 
disks, resulting in a larger total torque (see Fig. |9}. 

In Fig. [T7] the surface density for the a = 0.004 and a = 
0.008 simulations with the two different ortho-para gas ratios 
are displayed. 

The surface density is a little bit higher in front of the planet 
(4> > 180 and r < 1.0aj up ) and a little bit lower behind the 
planet (<p < 180 and r > 1.0aj up ) for 0:P(equi) compared with 



0:P(3: 1) for both of the viscosities shown, influencing the posi- 
tive corotation torque. This can also be seen in the F(r) distribu- 
tion (Fig. [Tol l, where the spike becomes larger, resulting in the 
observed higher total torque for the 0:P(equi) models. 

For the lower viscosity (a = 0.004), the density in the spiral 
waves exerted by the planet is higher compared to the high vis- 
cosity case (a = 0.008), which could also be observed for disks 
with a constant y (Fig. PHI , Also, the surface density near the 
planet is higher for the lower viscosity disks. These two features 
cause the higher torque in the a = 0.004 disks compared with 
a = 0.008 (for the same gas ortho-para ratios). 

When a = 0.004, the spiral waves in the 0:P(equi) disks are 
much denser compared with 0:P(3:1) disks, but it seems to be 
the other way around for a = 0.008. The difference in the spiral 
wave density influences the height of the Lindblad torques in 

Hg.na 

In total, the differences in disk structure between the two dif- 
ferent ortho-para mixtures are larger in a constant viscosity disk 
than in disks with an ff-viscosity. The reason for this may be in 
the complicated relation that arises when an a viscosity is used, 
because in this case v depends explicitly on the temperature. 



5. Zero-torque radius 

The torque acting on an embedded 20ME art h pl anets decreases 
with increasing distance from the central star (Bit sch & Klevl 
1201 lh . At some point in the disk, the torque becomes zero 
(zero-torque radius Rq) and the planet ceases to migrate due to 
planet-disk interactions. The zero-torque radius Rq is strongly 
dependent on planetary mass and on gradients in the disk. In 
addition the location of Rq depends on the mass of the disk 
(lBitsch& Klevl [201lh . and based on the results presented in 
Section 21 we also expect Rq to depend on the viscosity and the 
equation of state of the gas. 

In Figure [18] we show the torque acting on 20ME a nh plan- 
ets placed at various radii in fully radiative disks with different 
viscosities. The figure also features simulations with 0:P(equi), 
0:P(3:1), and a constant adiabatic index y = 1.4. The data of 
the si mulations with consta nt y and constant viscosity are taken 
from lBitsch & Kle"vl d201 lh . Because a high viscosity is needed 
to keep torques that are acting on a planet from saturating, high- 
viscosity disks should have their Rq at larger orbital distances 
than low-viscosity disks. This is indeed the case in Figure [18] 
where Rq for the a = 0.004 simulation is closer to the star than 
in either the a = 0.008 or v = 10 15 cm 2 /s simulation. For these 
viscosities, Rq is usually largest for the constant adiabatic in- 
dex simulations and is smallest for the 0:P(3:1) simulations. 
The most realistic treatment for the equation of state is likely 
0:P(3:1), so the other approximations tend to overestimate the 
location of Rq. Table[2]summarizes our findings. 

To understand these results in more detail, consider Figure s[TJ 
and [10] The torque acting on an embedded 20ME art h planet de- 
creases with increasing adiabatic index. At low temperatures, 
the adiabatic index for 0:P(3:1) is constantly increasing with 
increasing temperature over the range of temperatures of inter- 
est. For 0:P(equi), the adiabatic index decreases until y « 1.3 
and then increases until y — 1.67 at T ~ 20^. As the temper- 
ature decreases with increasing distance in the disk, the result 
that 0:P(equi) simulations have a larger Rq than in the 0:P(3:1) 
simulations is consistent with results from the previous sections. 
Using the same reasoning, the constant adiabatic index sim- 
ulations should have larger Rq than the 0:P(3:1) simulations, 
and usually have larger Rq than in the 0:P(equi) case, although 
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Fig. 17. Displayed are the surface density maps for 20ME a nh 
planets on fixed circular orbits in fully radiative disks. The vis- 
cosities and gas mixtures in the disks are (from top to bottom): 
a = 0.004 with 0:P(3:1), a = 0.004 with 0:P(equi), a = 0.008 
with 0:P(3:1), and a = 0.008 with 0:P(equi). The snapshots are 
taken in the heating/cooling equilibrium state of the disk. 
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Fig. 18. Torque acting on 20ME a nh planets in disks with different 
viscosities. Each plot features the two gas mixtures [0:P(3:1) 
and 0:P(equi)] and a constant y = 1 .4. From top to bottom the 
viscosity is: a = 0.004, a = 0.008, and v = 10 15 cm 2 /s. The 
simulations wi t h y - 1 .4 and a constant viscosity are taken from 
iBitsch&KievlfcOllI) . 



the latter is complicated by the break in monotonicity of the 
0:P(equi) y temperature dependence. 
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6. Summary and Conclusions 

In this paper we have studied the disk structure due to different 
viscosities (constant a and constant v, respectively) and chang- 
ing adiabatic indices. We then investigated how an embedded 
20MEarth planet migrates in such disks. 

The viscosity determines the structure of the disk. We find 
that the midplane density (temperature) decreases (increases) for 
increasing a-viscosity. The aspect ratio of the disk and the slope 
of the temperature profile f5j also increase for increasing viscos- 
ity. The simulation with a constant viscosity of v = 10 15 cm 2 /s 
seems to lie between the a = 0.005 and a = 0.006 simulations, 
which has allowed us t o mak e co mparisons with the ear lier sim- 
ulations of iKlev et alJ d2009t) and lBitsch & Kiev! l201Ch . 

The aspect ratio, Hj r, of the disk is approximately indepen- 
dent of the index y. Only for very small values, y<1.2,H/r drops 
due to the onset of convection in the disk. At very low viscosi- 
ties, the aspect ratio profile is roughly constant with small fluc- 
tuations. However, for higher viscosities, the aspect ratio profile 
has a complex shape, with a slight increase for increasing r at 
small disk radii, followed by decreasing ratio at large disk radii. 
A decreasing H/r profile with increasing distance is the result of 
the temperature dependence of the opacity. 

A change in the aspect ratio influences the migration of em- 
bedded planets in radiatively cooling disks. Although we find 
that the direction of migration can be inward or outward, the ef- 
fect of the H/r profile is broadly consistent with results from 
isothermal disks, where embedded pla nets migrate inward faster 
for smaller H/r dTanaka et al.l l2002t) . The slope of the tem- 
perature /3j also influences the torque acting on an embedded 
planet. In many theoreti cal formulae dMasset & Casolil |2010t 
iPaardekooper et ail 1201 ll) the parameter /3-r is used to calculate 
the Lindblad torque and the entropy-related corotation torque. 
A change in these parameters therefore changes the torque act- 
ing on an embedded planet. We note that in all cases, isothermal 
or radiative, vis cosity is requ ired to prevent corotation torques 
from saturating (Masset feOOll) . In fully radiative disks, radiative 
transpo rt/cooling acts to prevent saturation of the entropy related 
torque (IKlev et al.ll20Q9b . 

For a < 0.0025, we find that the torque acting on a 20ME art h 
planet is negative (inward migrat ion), consistent with planets 
embedded in an isothermal disk (Bit sch & Klev|[20Toh . In the 
low viscosity case, a low-mass planet seems to open a partial 
gap in the disk, such that corotation torques are weakened. For 
increasing viscosities (a > 0.003), the torque becomes positive 
and seems to stall for a > 0.005. Although the density, tem- 
perature, and aspect ratio profiles for the a = 0.005 simulation 
match quite well with the constant viscosity simulation of an 
unperturbed disk, there is a difference of about 50% in the to- 
tal torque acting on an embedded 20ME ar th planet, emphasizing 
further that the treatment for viscosity affects planet migration. 
Because the a-viscosity is dependent on the temperature in the 
disk (c s oc Vr), and because the presence of the planet modi- 
fies the temperature in its neighborhood, the torques will change. 



Overall, this results in a smaller torque acting on the embedded 
planet. 

In addition to viscosity, the adiabatic index influences the 
torque on the planet. In our test simulations with a constant adi- 
abatic index, the torque on a planet is about a factor of three 
larger for an equation of state with a fixed y = 7/5 than for a 
fixed 7 = 5/3. For the simulations in which the adiabatic index 
is allowed to vary with temperature, the torque acting on embed- 
ded planets in a disk with 0:P(equi) is larger than in simulations 
with 0:P(3: 1) for the same viscosity. This picture is complicated 
though by the connection between the viscosity and the adia- 
batic index. Higher viscosity results in higher midplane temper- 
atures, which affects the temperature-dependent adiabatic index. 
The adiabatic index in turn alters the structure of the disk by 
changing for example the H/r profile. As a result, the location 
of the zero-torque radius is moved in the disk. 

If viscosity and disk temperature is positively correlated 
with disk mass, we would then expect massive disks to sup- 
port outward migration of embedded planets. This may be a 
way of transporting cores to the outer nebula and maintain- 
ing low eccentricities, wh ich may be necessary to explain sys- 
tems such as Fomalhaut dKalas et al.ll2008t iBolev et alj|2012h 
and HR8799 ( Marois et al. 1 120081 12010ir However, the embed- 
ded planet would need to be prevented from undergoing rapid 
gas accretion during the migration. We suggest that this might 
be an opacity effect that is exacerbated by the high disk tem- 
peratures around A stars. The higher temperature will place the 
zero-torque radius at larger semi-major axes, and the naturally 
higher mass of A star disks will increase the dust opacity, for 
constant metallicity. This could prevent cooling and rapid gas 
capture onto the planetary core until much lower surface densi- 
ties are reached. 

In addition to affecting massive planetary cores, the zero- 
torque radius, could in principle serve as a location for collect- 
ing embryos, perha ps aiding in the formation of a large core 
dSandoret al.ll20lTt) . However, variations in viscosity as the disk 
evolves will cause the location of the zero-torque radius to also 
evolve. This could result in the loss of embryos and changes that 
result to inward migration, or could open the possibility of pro- 
viding new feeding zones for young planets. The overall influ- 
ence of the zero-torque radius on the formation of large cores is 
further complicated by collisions between planetesimals in the 
envisaged planet trap, which may influence the opacity. Such a 
change will affect the H/r profile, which ultimately influences 
the zero-torque radius. 

The 0:P(3:1) equation of state is the most realistic for disks 
explored here. This treatment for the equation of state also leads 
to the smallest zero-torque radii compared with both the con- 
stant y = 1.4 and the 0:P(equi) simulations. The effect is most 
pronounced in high-viscosity disks, and the inclusion of a proper 
equation of state is crucial for determining the location of possi- 
ble planet traps, where planetary embryos or cores can collect. 





0:P(equi) 


0:P(3:1) 


const, y = 1 .4 


a = 0.004 


1.4 


1.7 


1.75 


a = 0.008 


2.7 


2.8 


>3.5 


const, v 


2.25 


1.7 


2.45 



Table 2. Zero-torque radii (in aj up ) for planets in fully radia- 
tive disks with different viscosities and different thermodynamic 
treatments of the gas. 
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